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(Dated:) 

We propose a simple algorithm which produces high dimensional Apollonian networks with both 
small-world and scale-free characteristics. We derive analytical expressions for the degree distri- 
bution, the clustering coefficient and the diameter of the networks, which are determined by their 
dimension. 
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I. INTRODUCTION 

Since the ground-breaking papers by Watts and Stro- 
gatz on small-world networks 1] and Barabasi and Albert 
on scale- free networks 0, the research interest on com- 
plex networks as an interdisciplinary subject has soared 
[J Complex networks describe many systems in 
nature and society, most of which share three apparent 
features: power-law degree distribution, small average 
path length (APL) and high clustering coefficient. 

While many models [J have been proposed to 
describe real-life networks, most of them are stochastic. 
However, new deterministic models with fixed degree dis- 
tributions constructed by recursive methods have been 
recently introduced 0, 0, H H [H El Determin- 
istic models have the strong advantage that it is often 
possible to compute analytically their properties, which 
may be compared with experimental data from real and 
simulated networks. Deterministic networks can be cre- 
ated by various techniques: modification of some reg- 
ular graphs E3|j addition and product of graphs |l4j . 
and other mathematical methods as in |15|. Concern- 
ing the problem of Apollonian packing, two groups inde- 
pendently introduced the Apollonian networks |lfjl Il7j 
which have interesting properties like being scale-free, 
Euclidean, matching, space-filling and can be applied to 
porous media, polydisperse packings, road networks or 
electrical supply systems [l6| and may also help to ex- 
plain the properties of energy landscapes and the associ- 
ated scale- free network of connected minima 17]. 

In this paper we present a simple iterative algorithm to 
generate high-dimensional Apollonian networks based on 
a similar idea as that of the recursive graphs proposed in 
|12| . The introduced algorithm can concretize the prob- 
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lems of abstract high-dimensional Apollonian packings. 
Using the algorithm we determine relevant characteristics 
of high-dimensional Apollonian networks: the degree dis- 
tribution, clustering coefficient and diameter, all of which 
depend on the dimension of Apollonian packings. 

It should be pointed out that the concept of high- 
dimensional Apollonian networks was already introduced 
in E3 and ^3- I n these works, however, the emphasis 
is placed on two-dimensional Apollonian networks and 
their aim is to address the behavior of dynamical pro- 
cesses or provide a model to help understand the 
energy landscape networks 0, 0, UJ. Here, we focus 
on the producing algorithm, based on which we provide 
a detailed calculation of the topology characterization of 
high-dimensional Apollonian networks and we show that 
it depends on the dimension. 



II. THE CONSTRUCTION OF HIGH 
DIMENSIONAL APOLLONIAN NETWORKS 

From the problem of Apollonian packing, a two- 
dimensional example of which is shown in Fig. 1(a), An- 
drade et al. introduced Apollonian networks [16j which 
were independently proposed by Doye and Massen in |17| . 
Apollonian packing dates back to Apollonius of Perga 
who lived around 200 BC. The classic two-dimensional 
Apollonian packing is constructed by starting with three 
mutually touching circles, whose interstice is a curvilinear 
triangle to be filled. Then a circle is inscribed, touching 
all the sides of this curvilinear triangle. We call this the 
first iteration t = 1 and the initial configuration is de- 
noted by t = 0. For subsequent iterations we indefinitely 
repeat the process for all the newly generated curvilin- 
ear triangles. In the limit of infinite iterations, the well- 
known two-dimensional Apollonian packing is obtained. 

From the two-dimensional Apollonian packing, one can 
straightforwardly define a two-dimensional Apollonian 
network [ltl Il7| . where the vertices are associated to 
the circles and two vertices are connected if the corre- 
sWriiCo circles are tangent. Figure 1(b) shows how the 
network is translated from the two-dimensional Apollo- 
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(a) 




FIG. 1: (a) A two-dimensional Apollonian packing of disks, 
(b) Construction of two-dimensional Apollonian networks, 
showing the first iterations steps. Vertices introduced in the 
initial, first, second and third iterations have symbols Q, •, 
■ and □ respectively. 



nian packing construction. The two-dimensional Apollo- 
nian network can be generalized to high-dimensions (d- 
dimensional, d > 2) [TtJ associated with other self-similar 
packings [2(J]). Next a comprehensive account will be 
given. 

First, we address d-dimensional Apollonian packings 
which are constructed iteratively in a similar way as 
shown in Fig. 1(a). Initially, we have d + 1 mutu- 
ally touching d-dimensional hyperspheres with a curvi- 
linear d-dimensional simplex (d-simplex) as their inter- 
stice. In the first iteration one e?-hypersphere is added 
to fill the interstice in the initial configuration, such that 
the added d-hypersphere touches each of the d + 1 d- 
hyperspheres. The process is continued for all the newly 
created curvilinear d-dimensional simplex in the succes- 
sive iterations. In the limit of infinite iterations, the re- 
sulting pictures are d-dimensional Apollonian packings. 
If each d-hypersphere corresponds to a vertex and ver- 
tices are connected by an edge if the corresponding d- 
hyperspheres are in contact, then one gets (i-dimensional 
Apollonian networks. 



III. THE ITERATIVE ALGORITHM OF HIGH 
DIMENSIONAL APOLLONIAN NETWORKS 

In the iterative process for the construction of high- 
dimensional Apollonian networks at each iteration, for 
each new hypersphere added, d + 1 new interstices are 
created in the associated Apollonian packing which will 
be filled in the next iteration, when building networks, 
we can say in equivalent words that for each new ver- 
tex added, d + 1 new <i-simplices are generated in the 
network, into which vertices will be inserted in the next 
iteration. According to this process, we can introduce 
a general iterative algorithm to create high-dimensional 
Apollonian networks which is similar to the process that 
allow the construction of the recursive graphs introduced 
in [12. 

Before introducing the algorithm we give the following 
definitions. A complete graph Kd (also referred in the 
literature as d-clique; see |21|) is a graph with d vertices, 
where there is no loop or multiple edge and every vertex 
is joined to every other by an edge. Generally speaking, 
two graphs are said to be isomorphic if the vertices 
and edges of one graph match up with vertices and edges 
of other, and the edge matching be consistent with the 
vertex matching. 

We denote the (i-dimensional Apollonian network af- 
ter t iterations by A(d,t), d > 2, t > 0. Then the d- 
dimensional Apollonian network at step t is constructed 
as follows: For t = 0, A(d, 0) is the complete graph 
Kd+\ (or (d + l)-clique), and A(d, 0) has d + 1 vertices 
and edges. For t > 1, A(d,t) is obtained from 

A{d, t — 1) by adding for each of its existing subgraphs 
isomorphic to a (d+ l)-clique and created at step t — 1 a 
new vertex and joining it to all the vertices of this sub- 
graph (see Fig. 1(b) for the case d = 2). Then, at t = 1, 
we add one new vertex and d+ 1 new edges to the graph, 
creating d+ 1 new Kd+i cliques and resulting in the com- 
plete graph with d + 2 vertices, denoted -FQ+2- At t — 2 
we add d + 1 new vertices, each of them connected to all 
the vertices of one of the d + 1 cliques Kd+i created at 
t = 1 introducing (d + l) 2 new edges, and so on. 

Note that the addition of each new vertex leads to d+ 1 
new (d + l)-cliques and d + 1 new edges. So the num- 
ber of vertices and edges newly born at step U become 
L V (U) = (1 + rf) t,_1 and L e (t t ) = (1 + d) ti , respectively. 
Therefore, similarly to many real-life networks such as 
the World Wide Web, the (i-dimensional Apollonian net- 
work is a growing network, whose number of vertices in- 
creases exponentially with time. 

Thus we can easily see that at step t, the Apollonian 
network A(d, t) has 



N t = (d + l) + J2 L v(U)= {d+l J - + d+l (1) 
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vertices and 



TABLE I: Distribution of vertices and their degrees for A(d, t) 
at step t. 



\E\t = 



d{d+l) 



"I> e (*i) = 



d(d + l) (d+l) t+1 -d-l Num. vert. 



Degree 



ti=\ 



edges 

The average degree k t is then 



2(d+l) t+1 +d 3 + d 2 -2d-2 
(d + 1)' + d 2 + d - 1 



(2) 



(3) 



For large t it is approximately 2(d + 1). We can see 
when t is large enough the resulting networks are sparse 
graphs as many real-world networks whose vertices have 
many fewer connections than is possible. 



IV. RELEVANT CHARACTERISTICS OF HIGH 
DIMENSIONAL APOLLONIAN NETWORKS 

Below we will find that the dimension d is a tunable 
parameter controlling all the relevant characteristics of 
the d-dimcnsional Apollonian network. 



A. Degree distribution 

When a new vertex i is added to the graph at step U, 
it has degree d + 1 and forms d + 1 new (d + l)-cliques. 
From the iterative algorithm, we can see that each new 
neighbor of i generated d new (d + l)-cliques with i as 
one vertex of them. In the next iteration, these id + 1)- 
cliques will introduce new vertices that are connected to 
the vertex i. Let fci(t) be the degree of i at step t and 
Afci(i) be the increment of ki(t) minus ki(t — 1). Then 
Aki(t) of vertex i evolves as 



Afci(i) = h{t) - h(t - 1) = dAk4t - 1) (4) 

d + 1 and 



combining the initial condition ki(ti) 
Afci(ti) = d + 1, we obtain 

Aki{t) = (d+l^- 1 *- 1 

and the degree of vertex i becomes 



(5) 



h{t) 



■1+ A M*r, 



--(d+1) 



t m =u+\ 



d*-*« - 1 

d-1 



(6) 

The distribution of all vertices and their degrees at step 
t is given in Table[IJ from which we can see that the degree 
spectrum of the graph is discrete and some values of the 
degree are absent. To relate the exponent of this discrete 
degree distribution to the standard 7 exponent as defined 
for continuous degree distribution, we use a cumulative 
distribution P cum {k) = J2k'>k N ( k ' \t)/ N t ~ ^ 1_7 - Here 



d+1 
1 

d+1 
(d+1) 2 

(d+1)*" 3 
(d+l)*" 2 
(d+1)*" 1 



+ 1) 

(d + l)(E =orf J + l) 

(d + i)(E,=od J + i) 

(d + l)(d + l + l) 
(d + l)(d + l) 

d+1 



k and k' are points of the discrete degree spectrum. The 
analytic computation details are given as follows. 
For a degree k 



(d+1) 



it— l 



1 



d-1 



1 



there are (d+ l) 1 ^ 1 vertices with this exact degree, all of 
which were introduced at step I. 

All vertices introduced at time I or earlier have this 
and a higher degree. So we have 



N(k',t) = (d+l)+J2 L v(s) = 



(d+1)' -1 



k'>k 



As the total number of vertices at step t is given in Eq. Q 
we have 



(<*+!)(' 



jt-i 



1 



1 !-T 



d-1 



(d+l)'-l 



1 



^±2^1 + d+l 

_ (d + iy + d(d + i)-i 

~ (d + 1)* + d(d + l)-l 
Therefore, for large t we obtain 

= (d + l)'~* 

and 

ln(d+l) 



(7) 



7 : 



1 



lnd 



(8) 



so that 2 < 7 < 2.58496. We note that this result was 
already obtained by Doye and Massen in [l7| . 

Also, notice that when t gets large, the maximal de- 
gree of a vertex roughly equals to d 



t-i 



N, 



lnd/ ln(d+l) 



N, 



1/(7-1) 



B. Clustering distribution 

The clustering coefficient [l| of a single vertex is the 
ratio of the total number of edges that actually exist be- 
tween all k its nearest neighbors and the potential num- 
ber of edges k(k — l)/2 between them. The clustering 
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coefficient of the whole network is the average of all indi- 
vidual vertices. We can derive analytical expressions for 
the clustering C(k) for any vertex with degree k. 

When a vertex is created it is connected to all the 
vertices of a (d + l)-clique whose vertices are completely 
interconnected. It follows that a vertex with degree k = 
d + 1 has a clustering coefficient of one because all the 
(d + l)d/2 possible links between its neighbors actually 
exist. In the following steps, if a newly created vertex 
is connected to the considered vertex whose d nearest 
neighbors must also be linked by the new one at the same 
step. Thus for a vertex v of degree k, the exact expression 
for its clustering coefficient is 



C(k) 



d(d+l) 
2 



+ d{k-d-l) _ 2d(k - ^±1) 
~~ k(k - 1) 



fc(fc-i) 

2 



(9) 



depending on degree k and dimension d. Using this re- 
sult, we can compute now the clustering of the graph at 
step t, it is C t = St/Nt, where 



S t = (d+l) 



2d(D 



d+l \ 

2 I 



D (D - 1) 



E 

q=l 



2d(D q 



d+l 

2 



L v (q) 



D g (D q 



(10) 



where D = ^ + d and D q = (d + 1) ( ^l" 1 + 1 
given by eq. (6) are the degrees of the vertices created at 
steps and q, respectively. One can easily prove that for 
t > 7 and for any d > 2 the following relation hold true. 



C t > 



3d - 2 
3d - 1 



(11) 



Therefore the clustering coefficient of high dimension 
Apollonian networks is very large. Similarly to the power 
exponent 7 of the degree distribution, the clustering is 
also tunable simply by changing the value of control pa- 
rameter d. From Eq. (11), one can see that the clustering 
coefficient increases with d and approaches a limit of 1 
when d gets large. In the special cases d — 2 and d = 3, 
C equals to constant asymptotic values 0.8284 and 
0.8852, respectively. 



Diameter 



(d + l)-clique, and since any (d + 2)-clique during those 
steps contains at least the vertex created at step 1 , which 
is from the initial clique Kd+2 or A(d,l) obtained after 
step 1, thus the diameter is 2. 

Steps [|] + 2 to d + 2. In any of those steps, some 
newly added vertices might not share a neighbor in the 
original clique Kd+2 obtained after step 1; however, any 
newly added vertex is connected to at least one vertex of 
the initial clique Kd+2- Thus, the diameter is equal to 3. 

Further steps. Clearly, at each step t > d + 3, the di- 
ameter always lies between a pair of vertices that have 
just been created at this step. We will call the newly cre- 
ated vertices "outer" vertices. At any step t > d + 3, we 
note that an outer vertex cannot be connected with two 
or more vertices that were created during the same step 
< t' < t — 1. Moreover, by construction no two vertices 
that were created during a given step are neighbors, thus 
they cannot be part of the same (d + 2)-clique. Thus, 
for any step t > d + 3, some outer vertices are connected 
with vertices that appeared at pairwise different steps. 
Thus, there exists an outer vertex vt created at step t, 
which is connected to vertices ViS, 1 < i < t — 1, where all 
the is are pairwise distinct. We conclude that vt is nec- 
essarily connected to a vertex that was created at a step 
to < t — d — 1. If we repeat this argument, then we ob- 
tain an upper bound on the distance from vt to the initial 
clique Kd+2- Let t = a(d + 1) + p, where 2 < p < d + 2. 
Then, we see that vt is at distance at most a + 1 from 
a vertex in Kd+2- Hence any two vertices vt and wt in 
A(d, t) lie at distance at most 2 (a + 1) + 1; however, de- 
pending on p, this distance can be reduced by 1, since 
when p < [~|] + 1, we know that two vertices created at 
step p share at least a neighbor in Kd+2- Thus, when 
2 < p < [| 1 + 1, diam(A(d,t)) < 2(a + 1), while when 
[|1 + 2 < p < d + 2, diam{A{d, t)) < 2(a + 1) + 1. One 
can see that these distance bounds can be reached by 
pairs of outer vertices created at step t. More precisely, 
those two vertices vt and wt share the property that they 
are connected to d vertices that appeared respectively at 
steps t — 1, t — 2, . . . t — d — 1. 

Hence, the diameter increases by 2 every d+l steps. 
More precisely, we have the following result, for any d > 2 
and t > 2 (when t = 1, the diameter is clearly equal to 
1): 

diam{A(d,t)) = 2(L^| J + 1) + f{d,t) 



The diameter of a network characterizes the maximum 
communication delay in the network and is defined as 
the longest shortest path between all pairs of vertices. 
In what follows, the notations \x~\ and [^J express the 
integers obtained by rounding x to the nearest integers 
towards infinity and minus infinity, respectively. Now we 
compute the diameter of A(d, t), denoted diam(A(d,t)) 
for d > 2 : 

Step 1. The diameter is 1. 

Steps 2 to + 1. In this case, the diameter is 2, 
since any new vertex is by construction connected to a 



if t- L|=fJ(d+l) < [f] +1, and 1 



where f(d, i) 
otherwise. 

In the limit of large t, diam(A(d,t)) ~ ^j-, while 

N t ~ (d + l) t_1 , thus the diameter is small and scales 
logarithmically with the network size. 



V. CONCLUSION AND DISCUSSION 

In conclusion, we have introduced a general iterative 
algorithm to produce high dimensional Apollonian net- 
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works associated with high dimensional packings. The 
networks present the typical characteristics of real-life 
networks in nature and society as they are small-world 
and have a power-law degree distribution. We compute 
analytical expressions for the degree distribution, the 
clustering coefficient, and the diameter of the networks, 
all of which are determined by the dimension of the asso- 
ciated Apollonian packings. The high dimensional Apol- 
lonian networks introduced here, and the consideration of 
the method presented in Ref . |23 , allow the construction 
of high dimensional random Apollonian networks piij . In 
addition, it should be worth studying in detail physical 



models such as Ising models |24J and processes such as 
percolation, spreading, searching and diffusion that take 
place on the higher-dimensional Apollonian networks to 
know also their relation with the dimension. 
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